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Introduction 

The primary tool for three-dimensional transonic unsteady aerodynamic 
analysis is the XTRAN3S code. 1 In that code, finite difference methods are 
used to numerically solve the transonic small disturbance potential equation. 

Steady state far-field conditions are implemented on the boundaries of the 
computational region. 

Using steady far-field conditions causes disturbances that are incident 
on the boundaries to be reflected back into the computational domain. The 
reflected disturbances can cause significant errors in the calculated 
aerodynamic loading. This work presents characteristic boundary conditions 
that reduce the reflected disturbances. The boundary conditions are 
extensions of those developed for two-dimensional flows by Engquist and 
Majda2»3 a nd Whitlow. ^ Some representative results are presented to show 
the effectiveness of the characteristic conditions. 

Symbol s 

unsteady lift curve slope 
reference chord 

reduced frequency based on reference chord (o)C r /U») 
free stream Mach number 
time 

free stream speed 

streamwise, spanwise, and lateral coordinate directions 
wing root angle of attack 

Ms- /d f/s ^ 


c *a 

c r 

k 

Moo 

t 

Uoo 

x,y,z 



Y ratio of specific heats 

$ transonic small disturbance velocity potential 

<*> circular frequency 

Unsteady Transonic Small Disturbance Potential Equation 
The XTRAN3S code is used to solve the unsteady transonic small 
disturbance (TSD) potential equation 

+ A<(> x ) t = (E<J> X + F<{> x + G<jy) x 

+ (* y + H* x 4>y>y + (* z ) z (1) 

where «|> is a disturbance velocity potential normalized by U®c r , c r is a 
wing reference chord, and U» is free stream speed. The spatial coordinates, 
x, y, and z, are normalized by c r , and time, t, is normalized by or 1 , 
where <o is the circular frequency of unsteady motion. In (1), 

C = k 2 M 2 

oo 

A = 2kM 2 
00 

E = 1 - M 2 
00 

The constants F, G, and H may be defined either as proposed in ref. 5 

F - (y+X)m£ 

G = \ (y-3)M* 

H = - (y-I)H^ 

or as proprosed in ref. 6 
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F - - j [3-(2 -t)M^ 
G = - A M 2 

C « 

2 

H = - M ^ 

00 


In the original version of XTRAN3S, the boundary conditions on the outer 
edges of the computational domain are 


<f> = 

: 0, 

upstream 

♦x 

= o. 

downstream 

*y 

= o, 

wing root 

♦y 

= 0, 

side 

4>z 

= 0, 

top and bottom 


Formulation of the Characteristic Boundary Conditions 
To formulate the characteristic far-field boundary conditions, the cross 
flow terms are dropped from the TSD potential equation. That is, G = H = 0, 
and (1) reduces to 


c *tt + A *xt 


B<b +6 + <f> 

xx Y yy Y zz 


( 2 ) 


where B = E + 2F<J> X . 


Assuming B to be locally constant, the transformations 
x 


5 = 


7b 


A 2 

T - w x + ** 


2 R 

where D = 4C + , are used to transform (2) into the wave equation 


A = + <t> + A 

y Tt y yy zz 


A characteristic far-field boundary condition that satisfies (3) is 7 

d> + A • + ^ = 0 

Y x Y r r 

where r^ = $2 + y2 + z 2. in untransformed coordinates, (4) becomes 


(3) 

(4) 
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Allowing x to approach - » with y and z finite, the following first order 
plane wave condition at the upstream boundary is obtained 

+ 7B> *t " ♦x = 0 (6) 

Letting x -*• + « with y and z finite results in the downstream condition 

4 + 7b } +t + ^x ° 0 

At the top and bottom boundaries, z -*■ ± » with x and y finite, and the 
boundary conditions are 

| ± 4> z = 0 (8) 

where + and - represent the top and bottom boundaries, respectively. Outboard 
of the wing tip, y + + <*> with x and z finite, resulting in the condition 

§ * t + 4> y = 0 (9) 

A symmetry condition, as implemented in the original XTRAN3S code, is imposed 
at the wing root. The boundary conditions are summarized in fig. 1. 

Results 

• The boundary conditions were implemented in XTRAN3S and tested by 
calculating the unsteady force response for a flat plate rectangular wing with 
a pulse in wing root angle of attack a. All computations were made using the 
CDC CYBER 203 vector processor at the NASA Langley Research Center. The 
calculations were made for M« = 0.85 on a 60 x 40 x 20 grid in x,y,z that 
extended ±20c r in x, 6.66c r in y, and ±10c r in z. The extent of the 
grid in the z direction is less than that of the default grid, which extends 
±25c r in z.8 The wing tip was located at 3.33c r in the y direction. 



Using the pulse-transfer function technique of ref. 8, the frequency 
response function for the unsteady lift curve slope c* a was calculated 
with and without the characteristic boundary conditions. Fig. 2a shows a 
comparison of the lift response calculated using steady-state far-field 
conditions with that obtained using a kernel function method. 9 Below 
k « 0.6, the XTRAN3S results have spurious oscillations due to disturbances 
reflected from the boundaries. When the characteristic conditions, (6)-(9), 
were implemented, the reflected disturbances were small. Fig. 2b shows that 
the oscillations in the unsteady lift response have been eliminated by using 
(6)-(9) on the boundaries. 

The cause of the oscillations in the lift response obtained using steady 
state conditions can be seen in fig. 3. In the pulse-transfer function 
method, after a is increased to a maximum and returned to its starting value, 
calculation of the unsteady forces is continued until they return to their 
starting values. Fig. 3 shows that with steady state conditions implemented 
on the boundaries, the unsteady lift begins to decline toward zero but shows 
an increase beginning near a time when a disturbance has traveled 
approximately 20c r . This is the time required for a disturbance emanating 
from the wing to travel to the z boundaries and return. Thus, the increase in 
the lift at Ucot/c r = 20 and the spurious oscillations in the XTRAN3S lift 
response in fig. 2a are attributed to disturbances reflected from the z 
boundaries. When the characteristic conditions are used, the lift history 
shows a smooth decline toward zero after the initial transient. As a result, 
the lift response in fig. 2b has no spurious oscillations. 

For this study, no calculations were made for flows past wings with 
thickness. However previous calculations of two-dimensional transonic 


5 



flows^ indicate that the present characteristic boundary conditions would be 
effective in reducing reflected disturbances in those cases. 

Concluding Remarks 

Characteristic boundary conditions for the three-dimensional transonic 
unsteady small disturbance potential equation have been developed. They were 
Implemented in the NASA Langley version of XTRAN3S and tested for a flat plate 
rectangular wing with a pulse in angle. of attack. The unsteady lift response 
for calculations with and without the characteristic conditions shows that use 
of those conditions eliminates most of the disturbances reflected from the 
computational boundaries. 
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(a) Steady state boundary conditions 



(b) Characteristic boundary conditions 

Figure 2. Force response for a flat plate rectangular 
wing, = 0.85. 
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